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Abstract 

Results of our efforts to model decomposition in tundra are presented. 
Models were developed from field studies and specific laboratory analyses. 
Later, model predictions were evaluated against further field and laboratory 
data. Respirometric measures of decomposition and micro-organism meta- 
bolism are emphasized. 


Introduction 


The recognition that word models are an important initial step in the de- 
velopment of more rigorous mathematical models was explicit in the early 
stages of the International Biological Programme (Coloumbe & Brown 
1970). The U.S. Tundra Biome’s Decomposition Group and members of 
the International Tundra Working Groups on Decomposition and Micro- 
biology initiated development of mathematical models in June 1972 (see 
Bunnell & Tait 1974 and Bunnell & Scoullar 1975 for a discussion of histori- 
cal development). The objectives of these modelling efforts were essentially 
threefold: 


(i) to explain seasonal and annual variations in decomposition rates; 

(ii) to evaluate the influences on decomposition of moisture, temperature 
and substrate chemistry; 

(iii) to incorporate explicitly microbiological data and relationships into a 
model which would facilitate intersite comparisons within tundra and 
allow extrapolation beyond the tundra ecosystem. 


It was apparent that studies in weight loss were inadequate to address 
these three objectives. 
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Measures of weight loss integrate many simultaneously occurring pheno- 
mena and a large amount of variability within and between years is expected. 
The initial attempts at synthesis across a wide variety of tundra sites (Heal 
& French 1974) documented this tremendous variability. In developing the 
simulation models discussed here the approach was to relate specific pro- 
cesses to broad integrative measures within a biologically sound and mathe- 
matically rigorous framework. Thus, laboratory measures of specific pro- 
cesses, such as respiration and substrate utilization (Flanagan & Veum 1974; 
Flanagan & Scarborough 1974) could be related to field measures of weight 
loss (Bunnell & Scoullar 1975; Bunnell et al. in press and Bunnell et al. in 
press). Our premise was that the explicit recognition of component processes 
would account for the great variability in integrative measures (Bunnell & 
Tait 1974). In developing the models we have recognized two integrative 
approaches to decomposition studies. One holds that the total ‘organic mat- 
ter-microbial’ complex respires measurable CO, in relation to climatic and 
state variables; the evolved CO, closely reflects the decomposition process 
and rate. The second approach states that the microbial biomass, annual 
productivity and associated respiration reflect the turnover rate of organic 
matter and when assessed in relation to annual primary productivity, inform 
us of the status of organic accumulation. Each broad approach can be stated 
as a complex hypothesis. Models incorporating either of these hypotheses 
can be tested for a variety of ecosystems. We have explored the hypotheses 
in relation to the tundra ecosystem particularly at Barrow, Alaska. 

It is not the purpose of this paper to describe in detail the models built 
by the decomposition group of the U.S. IBP Tundra Biome Programme 
or the International Tundra Biome’s Decomposition Working Group. The 
paper presents precis of two models, GRESP and DCOMP, which con- 
stitute a still evolving framework within which we are exploring decomposi- 
tion processes and the possibility of developing an understanding of their 
role in ecosystems. The models were presented by Bunnell & Tait (1974), 
evaluated by Bunnell eż al. (in press a and b) and examined within an ecosys- 
tem perspective by Bunnell & Scoullar (1975). The approach is strictly syne- 
cological and differs clearly from the approach of Marr et al. (1963). Wealso” 
have employed data on growth kinetics of individual micro-organisms in 
the model of Marr et al. (1963) to simulate organism turnover rates in th 
tundra phyllosphere. This more autoecological approach has been used 
several investigators in other ecosystems (Babiuk & Paul 1970; Gray. 
Williams 1971). i: 

The general paradigms involved in most earlier models of decompositic 
(e.g. Jenny et al. 1949; Olson 1963; Minderman 1968; Jones & Gore, r 
press) were reviewed and compared to modelling within the Tundra Biomi 
of IBP by Bunnell & Tait (1974). The models discussed in the present pape 
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differ from earlier works in that the variations in decomposition rates are 
over short periods (daily) and explicitly acknowledge the influences of 
changing moisture, temperature and substrate quality. Further, they expli- 
citly relate microbial processes such as respiration to broad patterns of 
weight loss. Bunnell & Tait (1974) discuss the merits of explicit versus 
implicit treatments of theory. 

However, as have most earlier studies, the models GRESP and DCOMP 
ignore microfloral biomass. Other workers (e.g. Gray & Williams 1971) have 
brought attention to the difficulties encountered in relating field measures 
of annual microbial biomass transients and chemostat assessments of effi- 
ciency of microbial productivity and maintenance demands to field measures 
of primary productivity. Annual microbial maintenance demands calculated 
from chemostat determination of maintenance requirements may preclude 
growth of field microbial populations (Gray & Williams 1971). However, 
microbial growth is known to occur in the sites in question. Explanations 
offered for this paradoxical situation include: (i) underestimated produc- 
tivity; (ii) conservative or inflated estimates of yield coefficient and mainten- 
ance demand respectively; (iii) no real accumulation of organic matter in 
the system, but a decrease; and (iv) inflated estimates of microbial biomass. 

Here we test the model of Marr et al. (1963) using fungal measures from 
Barrow. 


“Laboratory, field and modelling methods 


“Microbial ecologists must brave a wavering episteme whenever laboratory 
studies of processes are related to field studies of the results of processes. 
Epistemological questions loom especially large when the relations are to 
be expressed unambiguously in mathematical models (Bunnell 1973). Gross 
as most field data are, they are ‘real’, integrated measures of decomposer 
activities and the influences of rate regulators upon them. We assumed 
respirometric measures of microbial activity to be especially ‘real’ and 
selected moisture, temperature and substrate quality as the most influential 
regulators of both respiration and weight loss. Using Gilson respirometers 
we examined influences of moisture, temperature and substrate quality on 
microbial respiration in above-ground substrates (Flanagan & Veum 1974; 
Bunnell et al. in press). These measures were incorporated in the models 
GRESP; DCOMP; as described by Bunnell et al. (in press, a and b). 

Standing dead tissues were calculated in the manner described below. 
Before presenting the calculations for substrate availability it must be 
emphasized that the figures presented are subject to spatial and temporal 
variations within and between specific subsites of the Barrow intensive study 
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site (Bunnell et al. 1975c). The study area referred to here is the same subsite 
of the Barrow intensive site referred to in Flanagan & Veum (1974) and 
Flanagan & Scarborough (1974). This subsite is an almost pure stand of 
Eriophorum angustifolium with randomly dispersed small (25-100 cm?) 
patches of Carex aquatilis. 

Approximately 100g m~? of new leaf tissue was produced annually on 
the site. The average standing crop of standing dead leaf tissue over four 
seasons was also about 100g m~? and on average was composed of 60, 32 
and 11 g¢m_~? of leaves that were one, two and three years old respectively. 
Four-year-old standing dead leaves were rarely encountered. The between- 
year variation inamountof standing dead m-? was relatively small (go gm-?— 
118gm-?). 

In the absence of lemming influences, or production of new leaf 
materials, approximately 30% of the weight was removed as leachate from 
the standing dead zone by autumn (fall) rains and the freeze-thaw cycle 
leaving approximately 70 g m~? yr~! to become new standing dead tissue (all 
figures given are averages derived from four years of observation). Because 
standing dead vegetation did not accumulate from year to year we have 
assumed that on average the decomposition losses from standing dead 
approximate the rate of production, about 70g m~?. Further, there is a loss 
of 70—60=10gm_~ of standing dead vegetation during the first growing 
season. The effects of intensive lemming grazing greatly alter the relation- 
ships outlined above (Bunnell et al. 1975c). However, during the period 1972 
to 1974 inclusive lemming numbers were generally much lower than the 
pronounced cyclic highs noted by Pitelka (1973). 

Annual and instantaneous moisture and témperature dependent weight 
losses in terms of CO, loss in respiration were determined earlier for stand- 
ing dead leaves (Flanagan & Veum 1974; Bunnell et al. in press, a). In this 
study these values are compared with calculated respiration for fungal main- 
tenance and production in standing dead leaf remains. The approach allows 
an evaluation of the model of Marr et a/. (1963) for predicting organic turn- 
over rates ina relatively simple natural system using empirical data concern- 
ing specific decomposer micro-organisms and their respiration in a native 
substrate. 


Results and discussion 


A. THE CHEMOSTAT MODEL 


By examining growth, respiration, production efficiency and main n 
demands of individual micro-organisms in dead leaves and incorporating” 
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leaf tissue. The model of Marr et al. (1963) states: 


dx = ds 


dx. 
where dt 'Stantaneous growth rate 


a=maintenance demand (gg-'hr-'), 
x= biomass of micro-organisms, 
Yy=production efficiency or yield coefficient (gg-! substrate), 


ds A 
and dt 'PStantaneous change in substrate. 


If growth ceases then 


dx 
dt ° 
and 
x d 
ed (19.2) 


From growth and respiration measures of batch cultures (e.g. Fig. 19.14 
and 19.1b) the values of ‘a and ‘y’ were calculated as 032X 10-4 gg! hr-! 
and 0:35 ¢g-' substrate respectively. The average standing crop of fungi 
(x in equation (19.1) and (19.2)) determined by the Jones & Mollison (1948) 


ds_ 032% 10-3 gg! hr-' x 0-0045 g g~! x 100g m~? 


dt o35eg7 


=41 X 10-4 g hr-! m-2 


Substrate utilization for maintenance over the active period of 120 days 
is thus 4-1 x 10-4 ghr-' m~ x 24 hr x 120— 1'2gm~, which leaves 7o—1-2 
or 68-8 g m-? of substrate for growth. 

To estimate the number of generations we can employ equation (19.3): 


y(S+-Nx)=Nx (19.3) 


where y and x are as equation (19.1), 
S represents the amount of available substrate, and 
N represents the number of generations. 
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Figure 19.1. Respiration and growth of a common leaf litter fungus on Barrow tundra 
growing on a mineral salt medium containing a known glucose concentration. Vertical bars 
represent standard deviations. 


Using a value of 68-8 for S, from equation (19.3), we obtain 


n AS 
N=xq=y) | 


0'35 g g7! X 68-8 g m~? | 
045gm (I—0-35gg7") | 


=82:3 generations 


The computed number of generations corresponds to a renewal every 
35 hours during the active season or a gross production of 37g fungal tissue) 
m~? yr-!. 

However 37 g m~? production of fungal tissue would result in 100% loss 
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of standing dead leaves annually. The average annual loss in field leached 
standing dead is approximately 10gm~? as measured per unit length of 
marked dead leaves. Weight losses in litter below, measured in litter bags, 
is of the same magnitude. On this basis we have assumed that 8-8 gm~-? 
are available annually for microbial growth in standing dead tissue, allowing 
for 15 generations of the average standing crop m~? of microbes or a produc- 
tion of 4:7 gm~?. These predictions of microbial production are compatible 
with the availability of substrate and allow us to calculate a yield coefficient 
(sensu Heal & MacLean 1974) of y'= 4:7/8-8=0-54. Unlike the yield coef- 
ficients of equation (19.1), and (19.2), y’ allows the recycling of microbial 
tissue towards its own production. There appears to be no published esti- 
mates of fungal production in the vegetation canopy with which we can com- 
pare our results, It may be pertinent to note that our estimates of generation 
time do not fall within the range reported for Parinkina (1974) for bacteria 
in tundra soils. 

To a limited extent we can evaluate the veracity and integrity of our 
assumptions by independent studies on microbial respiration in standing 
dead vegetation and weight losses from litter bags placed in the standing 
dead canopy. During the growing season of the first year following death, 
standing dead vegetation is diminished by about 10gm~? or about 14% 
(10g m-?/7og m~?). Flanagan & Veum (1974) predicted weight losses due 
to microbial respiration of 10:1 — 1 1:8% and reported measured weight losses 
from litter bags in the standing dead canopy of 9-9-+ 18%. These observa- 
tions suggest that microbial respiration alone accounts for a significant pro- 

_ Portion of the total weight lost, and do not contradict the assumptions made 
in equations (19.1), (19.2) and (19.3). The high proportion of weight loss 
due to microbial respiration is similarly compatible with the initial leaching 
losses (30%), largely complete by the beginning of the active growing season, 
and the relatively low abundance of invertebrate saprovores feeding above- 
ground (Bunnell et al. 1975c). In summary, we feel that careful examination 
of individual micro-organisms can lead to understanding of rates of decom- 
position and organic matter turnover at the ‘macro’ level of ecosystems. 
Admittedly the present study was relatively simple because of the single 
substrate and predominance of a few organisms (Flanagan & Scarborough 
1974), but it does suggest that autoecological studies of micro-organisms 


can contribute to the understanding of rates often viewed at an ecosystem 
level. 


B. MATHEMATICAL SIMULATION MODELS 


Flanagan & Veum (1974) utilized empirical relationships to show that weight 
loss attributable to C loss in respiration was often close to weight losses 
measured in litter bags, particularly in material 1-3 years of age. They 
employed generalized moisture and temperature responses of tundra decom- 
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Figure 19.2A. Response of microbial respiration in Barrow litter at o°C, (A); 5°C, (O); 
and 10°C, (@), to changes in moisture content. 


Figure 19.2B. Response of Barrow litters (Q=Carex aquatilis and @=total litter) at 
120% H,O (dry wt) to temperature variation. In A and B lines are plotted by hand and 
points represent the mean of ten samples each of which was read ten times during the experi- 
ments. 

poser systems such as those illustrated in Fig. 19.2. Among the tundra myco- 
flora the Q ,, value (averaged over 10°C increments from — 10°C to 25°C) 
was 4. Respiration continued down to —5-0°C; was depressed by tempera- 
tures above 30°C and was greatest at 25°C. In general cardinal moisture 
levels were 20%, 200-500% and 1,000% respectively for minimal, optimal 
and maximal levels. Moisture and temperature responses were inter- 
dependent over a wide range of both parameters. Usually, increasing tem- 
perature shifted the optimal moisture range towards the Y axis as in Fig. 
19.2a. From such data we constructed respiration response surfaces for litter 
and standing dead vegetation of varying substrate quality (e.g. Fig. 19.3). 
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Figure 19.3. Response surface of microbial respiration versus temperature and moisture 
in total litter at Barrow. Each grid point represents an average of 6 to 30 replicates (Flana- 
gan & Veum 1974). 


Having documented general patterns and developed empirical models 
the next step was to relate specified processes within a biologically sound 
and mathematically vigorous framework (Bunnell & Tait 1974). The first 
model formulated was given the mnemonic GRESP because it was based 
largely on experimentation with Gilson respirometry. This model assumes 
that the major influences on microbial respiration are moisture, oxygen, tem- 
perature, and substrate. It incorporates the influences of these factors on 
microbial respiration as deduced from laboratory data and integrates the 
effects over climatic variation within a specific substrate and season to pro- 
duce independent CO, based assessments of weight loss. These predicted 
losses generated from specific process studies can be related to the integrated 
weight loss as measured by litter bags. 


(i) The GRESP Model 


The rationale of the model is described in detail by Bunnell & Tait (1974) 
and Bunnell et a/. (in press, a). High and low moisture levels adversely in- 
fluence respiration rates. Equations (19.4) and (19.5) are first approximations 
of the moisture influence on microbial respiration as incorporated in the 
model. At moisture levels below optimum, respiration decreases according 
to equation (19.4), at high moisture levels respiration decreases according to 
equation (19.5). Thus, at low moisture levels the rate of moisture supply to 
respiring organisms is determined by coefficient a |. Coefficient a, determines 
the supply of oxygen to micro-organisms under conditions of high moisture. 
Both coefficients in equations (19.4) and (19.5) are assumed to follow 
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Michaelis-Menten relationships with the percentage moisture content of the 
decomposing substrates. 
M 


WL) FM (19.4) 


= water limitation reducing respiration, 

M=% water content on a dry weight basis, and 

a, = water content at which the substrate is half ‘saturated’ with 
water. 


M (19.5) 


where O,=oxygen availability, 
M=% water content on a dry weight basis, and 
a, water content at which 50%, of air spaces are water filled. 


The proportion of free channels is therefore 
aR: 
= CA +M a,+ M 


The substrate quality influences on respiration are expressed as 


a,» a," 10110 (196) 


where a, =the optimal respiration rate at 10°C when neither H,O nor O, 
are limiting, 
a,=Q,, coefficient, and 
T=temperature’C. 


The formal statement of the respiration model is: 


T—t10/10 


(19.7) 


where R(TM)=the respiration rate at temperature T and moisture level 
M. 


Estimating model coefficients 


Because of the model’s nonlinear form, estimation of coefficient values Gi 
through a,) requires nonlinear optimization techniques (Bunnell et al. in 
press a). Two determinants of the goodness of fit with respect to coefficients 
a, through a, were computed and designated r? and Thr The coefficient 
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of determination r? is not the Square of the correlation coefficient used in 
linear regression but it can be used in the same way since it determines the 
proportion of variation in the data accounted for by the model. The statistic 
Ta? determines the proportion of ‘accountable variation’ in the data ex- 
plained in the model. Its derivation is discussed by Bunnell et al. (in press, a). 
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Figure 19.4. Measured (a-d) and simulated (e-h) rates of microbial respiration in standing 
dead of Carex aquatilis (a and e=t year old; b and f=2 year old), and Eriophorum angusti- 
folium (cand g) and total litter (d and h) at Barrow, Alaska (from Bunnell ef al. in press, b). 


Table 19.1. ‘Goodness of Fit’ of the GRESP Function to Data from Different Substrates. (Modified from Bunnell er al. in press, a.) 


Number of 
Substrate type Observations 
Eriophorum (new) 54 
Eriophorum (1 year old) 79 
Eriophorum (2 years old) 65 
Carex (2 years old) 45 
Dryas (2 years old) 38 
Carex (1 year old) 96 
Total Barrow litter 346 


‘best’ parameters 
| 
(H,0%, dw) 


416 
74 
19 

213 

160 
59 

116 


a, 


(H,0% dw) 


416 
12,709 


3x10” 


822 
160 
13,300 
2,820 


a; 


(HIO, g-' hr>!) 


777 
184 
107 
296 
158 
204 
232 


a, 


(Qio) 


879 
2:56 
2:79 
218 
t75 
3:34 
374 


* An optimum solution had not been found. 


t The age of this material was in doubt and the number of treatment replicates at each observation point was only two 


Decomposition models 440 


In essence it approximates analysis of variance techniques in a nonlinear 
system, and permits sampling error to be removed from the fit between 
model and data. Measures of the goodness of fit of the model to data from 
different substrates are shown in Table 19.1. Examples of respiration re- 
sponse surfaces generated by the model are compared with measured sur- 
faces in Fig. 19.4. 

From Table 19.1, it can be seen: (i) that the model accounts for 71- 
96% of the accountable variation in respiration data; (ii) the model fits the 
respiration patterns of a number of diverse substrates equally well; (iii) that 
the coefficients a, through a, change between substrates and with age. In 
particular, it can be seen that as Eriophorum material ages, a, values decrease 
while a, values increase. As diminution and physical changes occur in the 
substrate the range over which moisture is not limiting to respiration in- 
creases. Bunnell & Tait (1974) predicted this phenomenon stating ‘the more 


Table 19.2. Results of sensitivity analyses on the GRESP model for total litter. 


Coefficient values that reduce the Maximum departure from 
goodness of fit by 1% of the best value best value 
Analysis number best Minimum Maximum Max-Min 
Coefficient 1 2 3 4 5 fit as % best as % best as °, best 
a, (H,0°a dw) $2 214 7# m2 m3 +84% 116% 
a, (H,O dw) 4.293 085 što? 108 O51 + 84% 161%, 
a, (WlO,/g/h) 193 345 rs 344380 + 68%, Ryn 
a, (Qio) 38 36 38 38 yi 3o +u% 14% 


physically reduced and pitted the substrate the smaller will be the “‘field 
capacity” and the larger the “retentive capacity” ’; (iv)a,, a measure of sub- 
strate quality, decreases with age of Eriophorum substrate. There is a near 
linear relationship between the best fit value of a, and the percent ethanol 
soluble fraction as determined by Flanagan & Veum (1974) and (v) a, is 
usually higher for younger substrates (8-79 for freshly dead Eriophorum) than 
for older ones (2°56 for one-year-old Eriophorum). 

Sensitivity analyses employed in evaluating the GRESP model are de- 
scribed by Bunnell et al. (in pressa). The results of sensitivity analysis (Table 
19.2) indicate that microbial respiration in litter is much more sensitive to 
substrate and temperature variables than to moisture levels. Estimates of 
the sensitivity of the respiration rate to changes in the coefficients a, through 
a, decrease in the order a,,a,,a, anda,, i.e. high moisture causing decreased 
O, availability has least influence on predicted respiration rates while tem- 
perature has the greatest influence. 

These observations are companionable with the theoretical construction 
of the model and, together with the general goodness of fit of the model to 
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data, indicate that the conceptual framework expressed by equation (19.7) 
is reasonably correct. 

It is important to emphasize that the model is presently limited in several 
ways, particularly in its applicability to field sites where temperatures are 
often above 10°C. Further laboratory studies will increase its temperature 
range. The model does not consider anaerobic decomposition processes, 
does not account for influences of leaching and microfauna, and in general 
may underestimate total weight losses occurring during decomposition. The 
model’s use is primarily one of estimating the proportion of weight loss attri- 
butable to loss of carbon via microbial respiration. For a specific substrate 
the closeness of this estimate to field measures of weight loss will reflect 
to some extent the magnitude of unaccounted microfaunal and leaching 
influences on decomposition. The model’s prediction is substrate 
specific. Influences of variation in substrate quality on rates of de- 
composition are addressed by the first extension of GRESP, the model 
DCOMP. 

It is a basic assumption of GRESP that the form of the respiration re- 
sponse is specific for a particular substrate. Acknowledging this assumption, 
the model may be employed to predict the respiration rate under any set 
of climatic circumstances. By employing the coefficients for Barrow total 
litter (Table 19.1), the temperature data for Barrow during 1971 (Weller 
et al. 1971) and moisture data of Flanagan & Veum (1974), a prediction of 
weight loss in total litter was made for the period 2 June through to 2 Sep- 
tember 1971. The model predicted a weight loss of 13-4%. The prediction 
falls midway between independent measures in litter bag weight loss experi- 
ments of 9-9%+1°8% (Flanagan & Veum 1974) and 15% (Benoit et al. 
1972). On the basis of this performance, it seems that the model can make 
realistic predictions for decomposition of a specific substrate during a season 
of known climatic conditions. 

In this way interannual differences in weight loss can be related to cli- 
matic variation explicitly. Similarly the model allows examination of in- 
traannual weight loss transients and inter- and intrasite variations. The 
influence of diurnal moisture and temperature patterns can be explored by 
the model. Bunnell & Tait (1974) have shown that there is a 12% increase 
in C loss from litter when computed using diurnal temperatures at two- 
hourly intervals as compared with average daily temperatures. Differences 
in diurnal climatic patterns may account for some between site variation 
in weight losses observed by Heal & French (1974). As expressed in equation 
(19.7) the model allows us to predict weight losses from a specific substrate 
(e.g. one-year-old litter or cellulose) placed in sites of different climate. Thus, 
we should be able to predict weight loss in Barrow litter placed, for example, 
in a tropical grassland. We could not predict weight losses in tropical grass 
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litter unless we constructed respiration temperature/moisture surfaces for 
that litter. 

Where substrate quality transients are steep during decomposition, due 
to fractions being removed by water flows, microfaunal activities, or to ‘soft’ 
substrate decomposition before recalcitrant substrate removal, then the 
GRESP model cannot be used to predict weight losses. 


(i) The DCOMP Model 


A second model DCOMP has been designed to assess decomposition rates 
under the influence of changing substrate quality. The DCOMP model 
assumes that a decomposing substrate contains a number of basic chemical 
constituents each having a specific temperature—moisture—chemically deter- 
mined respiration rate. Basically this concept is an extention of Minderman’s 
(1968) hypothesis, first to include temperature and moisture influences and 
second, to relate a microbial process, respiration, to weight loss. We assume 
that the total respiration rate of a substrate containing several chemical 
groups will be the sum of the respiration of microbes using the different 
groups. Bunnell et al. (in press, b) have discussed several formulations relat- 
ing substrate chemistry to microbial respiration and weight loss. The 
assumptions of the version discussed here are as follows: (i) each chemical 
category will have a temperature—moisture-respiration response surface of 
specific amplitude; (ii) each chemical response surface will have a similar 
shapeand (iii) any two surfaces will differ in amplitude by a constant chemi- 
cal-specific factor. The basic response surface used in DCOMP is the best 
fit of the GRESP model to Barrow total litter. It could just as well be the 
GRESP fit to cellulose placed in the field at Barrow. It should differ from 
a specific substrate’s temperature—moisture—respiration response surface by 
a specific and computable constant. 

The regression models used in calculating the respiration rate for specific 
chemical categories over varying temperature and moisture regimes are: 


RR=S rx; (19.8) 


i=ı 
where RR=total respiration, 
r,=respiration rate of compound i, and 
x=amount of compound i present. 


To include the influence of temperature (T) and moisture (M) we write 


RR(T.M)=J'r(T,M)x, (19.9) 


i= 
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To remove or account for the influence of T and M we write 
RR(T,M)=} b;GRESP(T,M)x; (19.10) 


In equation (19.10) the independent variables T and M are expressed 
as a surface such as those in Fig. 19.4, b=the ratio between the GRESP 
response surface for total litter and the response surface for compound i 
(i.e. a chemical specific scaler for the GRESP surface) and GRESP (T,M)= 
the underlying shape of the substrate specific response surface to T and M. 
For a substrate containing two different chemical groups, i.e. a ‘soft and 
a recalcitrant fraction we write: 


RR(T,M)=b,GRESP(T,M)x, + b, GRESP(T,M)x, (19.11) 


where RR(T,M)=total predicted respiration in ICO, hr~', 

b =the ratio between the amplitude of the GRESP surface 
for total litter and the amplitude of the respiration re- 
sponse surface for the ethanol soluble compounds, 

b =the ratio between the amplitude of the GRESP response 
surface for total litter and the amplitude of the respira- 
tion response surface for the recalcitrant compounds, 

GRESP(T,M)=the response surface in pICO,g-'hr-' of respiration 
versus T&M for litter from Barrow (see equation 
(19.7)), 
x,=amount of ethanol soluble compounds i, and 
x =amount of non-ethanol soluble compounds j 


Regression analyses employing field data are summarized in Table 19.3. 
The regressions are of the form 


In (% compound remaining)=b,+-b, (age in years) (19.12) 


The standing dead vegetation used in the analyses are E. angustifolium 
and C. aquatilis. Degrees of freedom are small because only one measure 
of chemical composition was associated with several measures of weight loss 
(Bunnell eta/. in press, b). Nevertheless, slopes of all regressions depart signi- 
ficantly from zero («= 0-05) which supports Minderman’s (1968) hypothesis 
of exponential decay for specific chemical constituents. The different slopes 
or decay rates for cellulose placed in different portions of the system reveal 
the effects of different microclimatic regimes. 

Results of analyses employing equation (19.10) are discussed in some 
detail by Bunnell eż al. in press, b. Whereas equation (19.12) is simply a 
least squares fit to field measures, equation (19.10) posits a specific relation- 
ship between microbial respiration and abiotic variables, and sums the re- 
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sponse of several component substrates. When the two broad chemical 
groups (equation (19.11)) were grouped by temperature and moisture classes 
significant relations (4=0-05) between respiration rate and chemical com- 
position were observed in only one third of the classes. 

These findings suggest that there is no simple underlying linear relation- 
ship between respiration rate and substrate composition. Broad relationships 
are, however, apparent. Regression analysis using equation (19.11) with 716 
measured respiration rates and nine substrate types yielded coefficients of 
3°8 and 0-76 for the ethanol soluble and non-ethanol soluble terms respec- 
tively. The associated t values were 5gand 15-2. With 7 d.f. we are confident 
that slopes of the regression differ from zero; X=o0-001. 

The model or equation (19.11) estimates the relative respiration rates 
to be 5:1 for ethanol soluble relative to the more recalcitrant fraction. This 
ratio approximates that observed for decomposition rates estimated from 
the standing dead data (0°69/0:12=5-75:1; see Table 19.3). 


Table 19.3. Regression coefficients of In (percent remaining) versus age in years. In 
(percent remaining) = ba +b, (age). (From Bunnell et al. in press, b.) 


Degrees 

of 
Substrate Environment bo b, r freedom « 
Cellulose Standing dead 45 —O-015 og8 2 0-05 
Cellulose Litter layer 458 —0-036 —o-97 2 0°05 
Cellulose Soil 460* —o-01g4 —0-997 1 O05 
Ethanol soluble Litter layer +412 —o-69 —o-go 8 oor 
Non-ethanol soluble Litter layer 456 —o12 —o68 8 0-05 


* In (100% )=4-605. 


Assuming that both compound groups contain 45% carbon the respira- 
tion rates may be converted to weight losses by separately multiplying the 
two additive terms of equation (19.11) by 


1zugmC rg substrate 
224HICO, ~~ o-gsugc 
and adding the two products. 

Using field temperature and moisture data from Barrow the annual 
decomposition rates of the two fractions were computed to be 0-66 and 
0'13 gg" which is very close to the decomposition rates observed in the 
field (0-69 and 0-12 g g-'; Table 19.3). Table 19.4 shows measured and simu- 
lated values for moisture-temperature-chemical specific rates of weight loss 
from E. angustifolium for the field season of 1973. Because it simulates res- 
piration and thus decomposition rates for specific chemical components, the 
DCOMP model also predicts chemical composition (Table 19.4). 
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Table 19.4. Measured and simulated values for moisture-tempera- 
ture-chemical specific rates of weight loss from Eriophorum angusti- 
folium. Measured values use changes in identified age classes. Simulated 
values use abiotic data from 1973 to project loss rates (from Bunnell 
et al. in press, b) 


Measured Simulated 
o, o, 
o © 


% loss per year: 
ethanol soluble 49 48 
non-ethanol soluble i 12 


Chemical composition 

of the substrate after one year: 
ethanol soluble 
non-ethanol soluble 

Total weight loss 


10 
go 
3! 


35y 


The results summarized in Tables 19.3 and 19.4 document that there 
is reasonably close agreement between simulated weight losses and substrate 
composition and the measured values of these. Presently we can assume that 
the biological principls of decomposition are realistically formulated in the 
model structure. Also, it appears that chemical composition and, therefore, 
rate of decomposition loss for specific substrate components can be predicted 
providing temperature and moisture influences on microbial respiration are 
incorporated in the model. Consideration of the rates of decomposition 
losses for specific substrates need not be restricted to carbon compounds as 
we have done here. There are specific N and P requirements for the decom- 
position of ethanol soluble and the more recalcitrant substrates. Relation- 
ships between C losses and N and P losses formulated in a GRESP or 
DCOMP model would almost certainly lead to useful insights regarding 
mineral release from decomposing (respiring) tissues. Consideration of N 
and P losses accompanying respiration and the possibility of predicting N 
and P losses and content of substrates from respiration data form a rationale 
for continuing to develop GRESP and DCOMP as predictive tools for 
understanding the role and significance of decomposition in ecosystems. 


Summary 


GRESP and DCOMP have examined the climatic and substrate influences 
on decomposition rates and represent a more synecological approach than 
that of the chemostat model discussed earlier. Still, the range and rate of 
decomposition determinable or predictable using either method are very 
similar. By two dissimilar routes we have arrived at a similar assessment 
of decomposition rates for plant leaf tissue at Barrow, Alaska. The simulation 
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models allow exploration of the role of major influencing factors on decom- 
position while providing a means of examining temporal and spatial variation 
in the rates of decomposition. Additionally these models allow us to predict 
overall rates of organic matter decomposition which if assessed in relation 
to productivity measures provide information on organic accumulation 
rates. The models ignore, amongst other variables, differences in microbial 
growth potentials. Rather, they sum all microbial growth potentials and in- 
corporate the universal micro-organisms responding only to temperature, 
moisture and substrate variables. They thus represent a restatement of 
Winogradsky (1887) ‘Everything is everywhere—the environment selects’. 

Conversely, the model of Marr et al. (1963) examines microbial growth 
kinetics in relation to substrate availability ignoring the climatic and sub- 
strate quality influences on microbial growth potential. Whether we use one 
approach or the other we arrive at the same assessment of decomposition 
rates in the standing dead leaf remains. When we use both approaches we 
arrive at a paradoxical situation where we are unable to document what are 
the major determinants of decomposition rates—microbial population 
decomposition potentials or the influence upon these of environmental 
factors. It would appear that both constraints apply, particularly at the extre- 
mities of decomposition, but that the function decomposition proceeds 
virtually unregulated in rate across a wide spectrum of micro-organism and 
climatic variables, at least in tundra. 

We may draw support for this statement by considering the breadth of 
the respiration versus moisture curve particularly at low temperatures (Fig. 
19.2a). Also the influences of high and low temperatures are diminished by 
low and high moisture levels (Fig. 19.2a), and low and high moisture influ- 
ences are counteracted respectively by high and low temperature levels 
(Flanagan & Veum 1974). To these observations we may add that micro- 
organisms in the tundra fall into three basic categories, i.e. psychrophiles, 
cold tolerant mesophiles and thermo tolerant psychrophiles (Flanagan et al. 
in prep.), and that individual micro-organisms contain different enzymes 
functioning optimally on specific substrates throughout a temperature range 
from —8°C to + 35°C (Flanagan & Scarborough 1974). Acknowledging 
these observations it is clear that the range for micro-organism activity is 
indeed wide and that the influences of climatic and substrate variables upon 
it are masked because subtle enzymes activities grading to optima within 
and between micro-organisms occur throughout broad temperature and 
moisture ranges. The result of such activities and relationships produces 
a very broad ‘optimal’ mode for decomposition which lends itself well to 
prediction by relatively simple modelling techniques. 
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